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Abstract 

Through a simple and exact analytical derivation, we show that for a particle on a lattice, there 
is a one-to-one correspondence between the spectra in the presence of an attractive potential V 
and its repulsive counterpart —V. For a Hermitian potential, this result implies that the number of 
localized states is the same in both, attractive and repulsive, cases although these states occur above 
(below) the band-continnum for the repulsive (attractive) case. For a 'PT-symmetric potential that 
is odd under parity, our result implies that in the "PT-unbroken phase, the energy eigenvalues are 
symmetric around zero, and that the corresponding eigenfunctions are closely related to each other. 
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Introduction: The energy spectrum of a quantum particle in an attractive potential V(r), 
in general, consists of discrete eigenvalues for which the eigenfunctions are localized in real 
space, and continuum eigenvalues with non square-integrable eigenfunctions. The energy 
spectrum for the corresponding repulsive potential — V(r) has only continuum eigenval- 



ues 



0,3- 



This situation changes dramatically when the particle is confined to a lattice or, 
equivalently, is exposed to a periodic potential. Indeed, repulsively bound two-atom states 
have been explored in detail since their experimental discovery in optical lattices 0, \a I and 
continue to be a source of ongoing work j^] in the context of the Bose-Hubbard model g, Q|. 
We note that in the Bose-Hubbard model, the interaction between the two atoms is short- 
ranged and is tuned via the Feschback resonance 3J. However, to our knowledge, the 
properties of single-particle states localized in the vicinity of a generic repulsive potential 
(defined below) have not been studied. In another area, localized states in parity + time- 
reversal (VT) symmetric one- dimensional lattice models, too, have been explored in recent 
years. These explorations have focused on the PT-symmetry breaking in the presence of 



attractive (real) on-site potentials with random PT-symmetric complex parts 111 ]. 

In this note, through a simple but exact derivation, we show that for a single particle on 
a lattice, there is a one-to-one correspondence between its energy spectrum in the presence 
of an attractive potential and the repulsive counterpart, and that the corresponding eigen- 
functions have identical probability distributions. For "PT-symmetric potentials that are 
odd under parity (and hence time- reversal), we show that if the PT-symmetry is unbroken, 
the energy spectrum must be symmetric around zero. 

One-dimensional Model: Let us start with the Hamiltonian for a particle on a one- 
dimensional lattice with only nearest-neighbor hopping energy J > 0, 

flo = -^X)(ctc i+1 + c] +1 c i ) (1) 

i 

where c] and Cj are creation and annihilation operators at site i respectively. The external 
potential is given by V = J2j Vj c j c j- We define the potential to be attractive provided 
^2jVj < and repulsive if is positive. Let = J2j fa,j\j) be an eigenstate of the 

Hamiltonian H + = H + V with energy E a where \ j) denotes a single-particle state localized 
at site j. The coefficients f a j obey the recursion relation 

- J [fa,j+l + fa,j-l] + Vjfaj = E a fa,j- (2) 
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We now consider the staggered wavefunction \<f> a ) = J^. f a j(—iy\j) ■ Using Eq.(j2J) it is 
straightforward to show that the staggered wavefunction satisfies the following equation 

H \(f> a ) = + |0 a ). (3) 

Thus, it is an eigenfunction of the conjugate Hamiltonian H_ = Hq — V with eigenvalue 
—E a . When V = 0, the energy spectrum is given by e k = — 2J cos(ka) and represents 
the well-known continuum band from — 2J to 2J where a is the lattice spacing. In this 
trivial case, indeed the eigenfunction \ip k ) = J2j sin(fcj) |j) and its staggered counterpart 
\4>k) = Z~2j s i n K 71 ~ \j) have energies ±e k respectively. 

Our result shows that if an attractive external potential V has n bound states below 
its continuum with energies E m (m = 1, . . . , n), then the corresponding repulsive potential 
—V must have an equal number of bound states above its continuum with energies —E m . 
Since the staggered wavefunction \<p a ) varies over the lattice length-scale a, it is energetically 
expensive and ill-defined in the continuum limit a — > 0. Physically, in the continuum limit, 
the absence of lattice-site scattering centers makes it impossible for a particle to localize 
near the repulsive potential. However, on a lattice, the probability distributions for the two 
states - a localized bound state \if) a ) with energy E a < — 2J in an attractive potential and the 
localized bound state \<j) a ) with energy — E a > +2 J in the repulsive potential - are identical. 
As a concrete example, we numerically obtain the spectrum for a lattice with N = 29 sites 
and a quadratic potential that vanishes at the ends, V m = A(m — 1)(N — m)/A^Q, where 
m = 1, . . . , N, N = (N + l)/2 is the center of the lattice and V/v = A. Figure [1] shows the 
ground state wavefunction ipGm for the attractive case, A/J = —0.5, (left panel) along with 
the highest-energy state wavefunction <p m for the repulsive case, A/J = +0.5 (right panel). 
It is clear that the two wavefunctions are related by <p m = (— l)" 1 " 1 " 1 ^,™- 
Two-particle Case: We can generalize this result in a straightforward manner to treat in- 
terparticle interaction U = £V . Ui-jhihj where the on-site number operator is given by 
fi>i — c \ c i- 111 the two-particle sector, the recursion relation satisfied by the relative-coordinate 



wavefunction is given by 



BE 



- Jk [< m+1 + < m _ x ] + C/(r m )< m = El < m . (4) 

Here —it /a < K < n/a is the lattice momentum associated with the center-of-mass of the 
two particles, Jk = J cos(Ka) is the effective hopping energy, r m = am = a{i — j) is the 
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FIG. 1. (color online) (a) The left panel shows the dimensionless ground-state wavefunction ipGm 
for an attractive quadratic potential V m = A(m — 1)(N — m)/iVg where N = 29 = (2N$ + 1) is 
the lattice size and A/J = —0.5. As expected for a quadratic potential ground-state, V'G.m is a 
Gaussian with width xq = a(A r gt/|A|) 1/ ' 4 ~ 4.61. (b) The right panel shows the dimensionless 
highest-energy state wavefunction cf> m for its repulsive counterpart with A/J = +0.5. We see that 
the 6 m is indeed the staggered version of the ground-state wavefunction (f>G,m- 



distance between the two particles on the lattice located at sites i and j, and U{r m ) is 
the real-space interaction between the two particles. Note that for a non-local interparticle 
interaction U(r m ), multiple bound-state solutions are generic, although, in the context 



of the Bose-Hubbard model, only one [3| or two js] have been discussed. If ip^ is an 
eigenfunction of the Hamiltonian H + U with energy E^, Eq. H] implies that the staggered 
wavefunction 0^ denned by 4>a( r m) = ( — l) m ^( r m) is an eigenfunction of the conjugate 
Hamiltonian H — U with energy —E^. 

Two-particle bound states in the presence of on-site and nearest-neighbor repulsive 
density- density interactions on a lattice have been extensively investigated 

m. 

Our 

derivation shows that they are a generic feature of any density-density interaction on a lat- 
tice, and this result is true for square lattices in higher dimensions. Note that the quantum 
statistics of the particles only constrains the relative wavefunction ipa( r m) to be odd (spin- 
less fermions) or even (bosons or spin-singlet fermions) under parity; however, it does not 



affect the one-to-one correspondence between the spectra for the two Hamiltonians Hq ± U . 
Thus, two-atom bound-states with attractive and repulsive interactions in optical lattices 
(bosons) [3J, the donor and acceptor impurity levels in semiconductors (fermions) jgj], as 
well as the localized phonon modes (collective bosonic excitation) (9), [lo| around a soft or 
stiff impurity can all be thought of as manifestations of the correspondence between spectra 
for H + and 

VT Symmetric Potential: The mapping between the two Hamiltonians H + and H_ is valid 
independent of the properties of the potential V including its Hermiticity; the on-site po- 
tential elements Vj may be complex. However, for a PT-symmetric potential that is odd 
under parity (and hence, time reversal), V* = —Vj = V-j, it follows that = H- where 
* denotes complex conjugation. Therefore, it follows from H + \ip a ) = E a \ip a ) that the wave- 
function |^*) = fa,j\j) i s an eigenstate of the conjugate Hamiltonian H_ with eigenvalue 
In the continuum limit, it has been shown that a wide class of such potentials, in- 



cluding V(x) = ix 3 and V(x) = i sin 2n+1 (x) have purely real energy spectra 12|, Il3 |. If 
the PT-symmetry is unbroken, E* a = E a , then it follows that H_\cf) a ) = —E a \(p a ) and 

H.\r a ) = +E a m. 

This explicit construction of wavefunctions with equal and opposite energies implies that 
for any arbitrary "PT-symmetric potential that is odd under parity, if the VT symmetry 
is not broken, the energy spectrum must be symmetric around zero. It also shows that 
the corresponding wavefunctions in the two cases have components that are simply related: 
[+E a ,f*j] <H- [— E a , f a ,j{— As an example, we consider the simplest "finite lattice" 
with 2 points. (Our result is equally applicable to a finite lattice.) The Hamiltonian in this 
case is given by H- = —Jo~ x + i^o z where (o~ x ,o~ z ) are the Pauli matrices in the site- index 



space 14j and a real 7 ensures that the potential is odd under parity as well as time-reversal. 



The eigenvalues in this case are given by E± = ±a/J 2 — 7 2 . Thus the PT-symmetry in this 
case is not broken as long as 7 < J. The corresponding (unnormalized) eigenf unctions 



are given by 14\ 



where 6 = arctan(7/ a/ J 2 — 7 2 ) is real when 7 < J. Therefore, in the PT-unbroken phase, 
the eigenvectors for positive and negative energies indeed are related by f-j = (— 
where j = 0,1. 



Conclusion: Our result, through a one-to-one mapping between attractive and repulsive 
potentials on a lattice, shows that localized states in repulsive potentials are ubiquitous. 
These states can be explored via local measurements. In contrast to the bound-states with 
energies below the continuum band, these localized states with energies above the continuum 
band will decay into the continnum states. They may thus provide a useful spectroscopic 
tool in optical lattices as well as engineered electronic materials with a small bandwidth. 
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